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ABSTRACT 

The difference between the spectra of iron and other cosmic rays 
is interpreted in terms of two source mechanisms. One mechanism, 
possibly acceleration at neutron star surfaces, produces the iron and 
another is responsible for the rest of the primary nuclei. Within this 
model, high energy observations could determine whether secondary nuclei 
are produced in the sources or in the interstellar medium. 
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Recent experimental data on the composition of the nuclear cosmic . . 
radiation (l,2i3V4) inthe 1 to 100 Ge¥/nucleon region has raised 
questions on the previously accepted notions that-, all cosmic ray^ spectra 
are the same at high energies. These data indicate- that the- .ratios of 
both medium-to-iron nuclei and secondary-to-primary cosmic rays decrease 
with increasing energy i (Primary nuclei are directly accelerated in .. 
the cosmic ray sources and secondary nuclei are produced predominantly 
by spallation reactions of the primaries with matter between the sources, 
and earth). In this paper we wish to explore the implications on source 
composition and cosmic ray propagation of data obtained from a balloon- 
borne ionization spectrometer (5). We will show- that this data implies 
that iron nuclei have a~ different origin than the rest of the primary 

cosmic rays. ■ ■ ■ ' ' " . ■ ■ 

Spectra of protons, alpha particles, carbon, oxygen, 10 £ .Z s ,.14, 
and iron group nuclei are 'shown, in Figure 1. In limited energy ranges - 
thdse spedtra can be represented by power laws cp E ,: where 7 iS: the , ■ 
spectfS'I index'. The heavy solid and heavy dashed lines represent . 
maximum likelihood fits to -the -medium (C,N,0) and. Fe groups- in the. 

3 to 50 GeV/nucleon region. ' The differential spectral indices of these- 
groups are 2.64 + 0, 04 and • 2.12 |f Oi 13 respectively , (2) . The power . , 
law fit to the spectrum of medium (M) 'nuclei also fits the, 10 ^ ,Z ^ 14 
group; in addition the same power law fits the proton and alpha 

particle spectra in the 5 to 50 GeV/nucleon .fegibn.:, although at 

higher energies the spectra appear somewhat steeper. Thus, within the 
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experimental uncertainties, all primary nuclei with Z ^ 14 have similar 
spectra in the same energy per nucleon region. Note, however, the 
large difference of 0.52 + 0.15 between the spectral index of the iron 
group and that of primary nuclei lighter than iron. 

In Figure 2 the ratio of the secondary nuclei Li, Be and B (the 
light or L nuclei) to their immediate progenitors C, N and 0 is plotted 
as a function of energy per nucleon. The spectral index of the L 
group is 2.78 + 0.07 (2). The difference in spectral index between the L 
and M groups is 0.14 + 0.08 from about 1 to 20 GeV/nucleon in agreement with 
a similar difference reported by Smith et al. (4). Both 
measurements indicate that secondary nuclei have steeper spectra 

than their primary progenitors. However, the difference between the 
spectrum of iron and the spectrum of the other primary nuclei is sig- 
nificantly larger than the difference between primary and secondary 
nuclei. A similar conclusion may be obtained from iron and its 
secondaries (2) although with less compelling statistical significance. 

The difference in spectral index of primary and secondary nuclei 
can be best interpreted in terms of energy dependent propagation of 
cosmic rays- in the interstellar medium (1, 6). In a steady state 
model with exponential distribution of path length (e.g. Ref. 7), the 
ratio of the fluxes of L nuclei to M nuclei may be calculated from 

'Pm " 


1 + 'PlH ^ ^ML 'PFe 

^T.HT. 'Pm ^FeL '“’m J 


( 1 ) 
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where cp is a cosmic ray flux at earth, and the subscripts LH and Fe 
denote nuclei with (10 < Z ^ 14) and (25 < Z), respectively (L and 
M are defined above). The quantity X is the e-folding path length of 
the exponential distribution (or the mean escape path length of cosmic 
rays from the galaxy); is the nuclear destruction path length 

of L nuclei; and X^^^, X^^^^ and X^^^^ are the fragmentation path lengths 
of M, LH and Fe nuclei into L nuclei. 

-2 

Based on cross sections from Ref. 8, we have that X^^ = 13 g cm , 

= 19 g cm"^, Xlhl = 29 g cm"^ and Xp^L = 33 g cm"^. We take the 
observed flux ratios from Figures 1 and 2 and Ref. 2 as follows: 
cPl/cPm = 0.23E"°’^'^i°'°^, cplh/PM = 0*28 and cppe/'PM = 0.036E 0-52+0.15^ 
The path length X as determined by solving equation (1) is 

plotted as a function of energy/nucleon in the lower part of Figure 
3. At 1 GeV, X = 5 g cm"^, a value consistent with previous calcula- 
tions (7) of cosmic ray fragmentat±»n. The shaded area represents the 
uncertainty in X introduced by uncertainties in the spectral indices 
of both and the error bars represent uncertainties from 

CD /cD alone. The uncertainty in cp„ has almost no effect on X below 
M Fe M 

~ 10 GeV /nucleon and about a 10% effect at 40 GeV. 

Consider now the iron-to-medium ratio. In the steady state model, 
the source ratio (Fe/M)^ may be calculated from 


(Fe/M), - — 1 + 


( 2 ) 
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where Xj^ and Xp^ are the destruction path lengths of medium and iron 
nuclei and cp is the flux of medium nuclei of secondary origin. The ■ 

Mg 

ratio cpw /cp„ is energy dependent because of the energy dependence 

Mg M 

of the iron- to -medium ratio. However, because of the relatively small 
fragmentation of iron into the M group, we shall use the constant 
value cpj^ /cpj^ = 0.16 previously determined (7). 

The iron-to-medium source ratio from equation (2) is plotted in 
the upper part of Figure 3. As before, the shaded area represents the 
uncertainty in (Fe/M)g introduced by uncertainties in both 
cPf /cp^, whereas the error bars result from alone. 

The source ratio (Fe/M)g is about 0.08 at 1 GeV/nucleon, This 
value is in good agreement with ratios previously calculated from 
energy independent studies: (Fe/M)g = 0.11 (8), (Fe/M)g 0.1 (7), 
(Fe/M)g 0.09 (9). At higher energies, however, (Fe/M)^ exhibits a 
significant departure from a constant which cannot be explained by 
propagation effects alone. That propagation effects have only a small 
influence on (Fe/M)g can be seen from the fact that most of the uncer- 
tainty in (Fe/M)g comes from uncertainties in the ef-fect of 

CPl/c% is small. 

Straightforward evidence against a constant iron-to-medium source 
ratio also comes from the comparison of the observed energy dependent 
ratio cpp ^ 0 .036E^' with the low energy value of (Fe/M)^. 

Because more iron nuclei are broken up during propagation than medium 
nuclei, the iron-to-medium ratio at the source must always exceed the 
iron-to-medium ratio at earth. The latter, however, increases with 



5 


increasing energy and becomes greater than 0.1 somewhere between 5 
and 15 GeV/nucleon. This fact clearly requires an additional mechanism 
for iron production at higher energies. 

From these considerations we conclude that iron is produced by a 
different source mechanism than all other primary cosmic rays. A reasonable 
possib iLity would be the acceleration of iron nuclei in pulsars since 
the surface of the neutron stars in these objects are believed to 
consist principally of iron (10) . Because all primary nuclei except 
iron appear to have the same spectrum, they are interpreted as having 
a common origin. They could then be produced at any of the proposed 
sites of cosmic ray acceleration such as supernova envelopes (11) and 
supernova remnants (12). 

Let us now examine the consequences and predictions of such a 
two component model. The differences in spectral index between L and 
M nuclei is within errors the same as the reported difference (2) 
in index between iron and its fragmentation products. This result 
would imply that most of the fragmentation takes place in the inter- 
stellar medium and not in the sources. However, data on the energy 
dependence of the (17 ^ Z ^ 25)/(Fe + Ni) ratio by Webber et al. (4) 
appears to indicate that this ratio decreases with increasing energy 
above 2 GeV/nucleon much more rapidly than the L/M ratio 
shown in Figure 2. If this difference is upheld by future data, we 
shall be able to conclude that most of the cosmic ray fragmentation 
takes place in the sources and not in the interstellar medium. We 
should note, however, that in this case the cosmic rays should sample 
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a sufficiently low average density in the interstellar medium. This 
could be achieved if cosmic rays preferentially avoid interstellar 
clouds . 

Consider now in Figure 1 the proton spectrum above 50 GeV and 
the alpha particle spectrum above 10 GeV/nucleon. The maximum likeli- 
hood method yields spectral indices of 2.75 + 0,03 and 2.77 + 0.05 for 
protons and alphas, respectively (13). As discussed above, however, 
the medium spectrum with index 2,64 fits both the protons and alphas 
below 50 GeV/nucleon. Thus, there seems to be evidence for a steepening 
of the proton and alpha particle spectra, A consequence of our model 
is a similar steepening in the spectra of all nuclei with Z ^ 14, This 
spectral change beyond about 50 GeV/nucleon could either be due to 
propagation effects in the interstellar and interplanetary media, or 
it could be produced in the sources. Observations may differentiate 
between these two possibilities J the spectrum of iron should steepen 
above 50 GeV/nucleon if we are observing a propagation effect; it 
should remain the same as below 50 GeV/nucleon if the steepening is 
produced in the sources of the Z ^ 14 nuclei. It should be noted, 
however, that if iron is accelerated at the surfaces of neutron stars, 
its spectrum may have a different high-energy cutoff due to photonuclear 
disintegration (14), The cutoff energy, however, depends greatly on 
the pulsar model, and no firm predictions can be made. 

One final prediction of the present model concerns VVH nuclei (15). 
Since these nuclei are believed to be produced by neutron capture 
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processes and probably are not present on neutron star surfaces, they 
should not be generically related to iron nuclei. Their spectrum 
should therefore differ from the spectrum of iron. 

In summary, we find that the increase of the ratio of iron-to- 
medium nuclei with increasing energy cannot be explained by propagation 
effects alone, and it appears to require an additional source or 
acceleration mechanism for iron at high energies. Based on a model 
in which iron comes from a different source than the rest of primary 
cosmic rays, we have made predictions of several observable effects. 

It is hoped that future high energy cosmic ray experiments, such as 
those planned for NASA's High Energy Astronomical Observatory Satellite 
(17), can make these observations. 


R. Ramaty 

V. K. Balasubrahmanyan 
J. F. Ormes 

Laboratory for High Energy Astrophysics 
Goddard Space Flight Center 
Greenbelt, Maryland 20771 
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FIGURE CAPTIONS 

1. Intensities of primary cosmic rays from protons to iron. The data 

(16) 

are taken from the Ref. 5 and from Ormes and Webber 

2. The ratio of light to medium nuclei. The data are taken from the 
references indicated in the Figure. 

3. The source ratio of iron to medium nuclei, (Fe/M)g, and the escape 
path length of cosmic rays in the galaxy. The shaded areas and 
error bars and are discussed in the text. 
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